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Abstract Prion diseases, traditionally referred to as transmissible spongiform en- 
cephalopathies (TSE), are invariably fatal and highly infectious neurodegenerative 
diseases that affect a wide variety of mammalian species, manifesting as scrapie in 
sheep, bovine spongiform encephalopathy (BSE, or 'mad-cow' disease) in cattle, and 
Creutzfeldt- Jakob disease (CJD), Gerstmann-Strussler-Scheinker syndrome(GSS), 
Fatal familial insomnia (FFI) and Kulu in humans, etc. These neurodegenerative 
diseases are caused by the conversion from a soluble normal cellular prion pro- 
tein (PrP*') into insoluble abnormally folded infectious prions (PrP^*^). The hy- 
drophobic region PrP(109-136) controls the formation into diseased prions: the nor- 
mal PrP(l 13-120) AGAAAAGA palindrome is an inhibitor/blocker of prion diseases 
(Mol Cell Neurosci 15: 66-78), and the highly conserved glycine-xxx-glycine motif 
PrP(l 19-131) can inhibit the formation of infectious prion proteins in cells (7 Biol 
Chem 285: 20213-20223). This article gives detailed reviews on the PrP(109-136) 
region and presents the studies of its three-dimensional structures and structural 
dynamics. 
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1 Introduction 



Prion diseases such as Creutzfeldt- Jakob disease (CJD) in humans and bovine spongiform 
encephalopathy (BSE or 'mad-cow' disease) in cattle are invariably fatal neurodegenera- 
tive diseases. Prions differ from conventional infectious agents in being highly resistant 
to treatments that destroy the nucleic acids found in bacteria and viruses. The infectious 
prion is thought to be an abnormally folded isoform (PrP^'^) of a host protein known as 
the prion protein (PrP*^). The conversion of PrP*-^ to PrP^'^ occurs post-translationally and 
involves conformational change from a predominantly a-helical protein to one rich in 
/3-sheet amyloid fibrils. Much remains to be understood about how the normal cellular 
isoform of the prion protein PrP*^ undergoes structural changes to become the disease as- 
sociated amyloid fibril form PrP^'^. The hydrophobic domain of PrP*^( 109-136) is highly 
conserved, containing a palindrome and the repeats of the GxxxG protein-protein inter- 
action motif (two glycines separated by any three residues; please note that the minimum 
number of residues to form fibrils should be 5 [1]). It is reported that the palindrome 
AGAAAAGA is an inhibitor/blocker of prion diseases [ L, 2J and the glycine-xxx-glycine 
motif GAVVGGLGGYMLGis also an inhibitor of prion diseases Sims. The alterations 
of residues in AGAAAAGA and GAVVGGLGGYMLG will drastically affect the abihty 
of cells and lead to the amyloid fibril formations (e.g. A 117V will cause the Gerstmann- 
Straussler-Scheinker prion disease, and the numerous mutants in [3]). Our computational 
results also confirm the amyloid fibril formation ability of the PrP( 109-1 36) hydrophobic 
region (Figure 1). As shown in Figure 1, if energy is less than the threshold energy -26 
kcal/mol then amyloid fibril is formed in the corresponding region of residues, thus, the 
palindrome segment PrP(l 13-120) and the GLGGY segment PrP(124-128) can be firmly 
confirmed having a strong amyloid fibril formation property. This paper will give detailed 
reviews on the PrP(109-136) region from the 3D molecular structure (Figure 2) point of 
view and presents the studies of its molecular structural dynamics. The rest of this paper is 
arranged as follows. Section 2 will give a survey of the research works on AGAAAAGA 
and present the 3D structure of amyloid fibrils in the AGAAAAGA segment. Then, in 
Section 3, the analysis of the glycine-XXX-glycine motif GAVVGGLGGYMLG inhibit- 
ing prion diseases will be done from the molecular structural point of view. Section 4 will 
make some concluding remarks on the hydrophobic PrP(109-136) region. 

2 Studies on the PrP(113-120) AGAAAAGA 
2.1 A Survey of the Research Works on AGAAAAGA 

The highly conserved hydrophobic palindrome AGAAAAGA PrP( 113-120) has been 
considered essential to PrP conformational conversion. Firstly we give a survey of the re- 
search works on prion AGAAAAGA listed in the PubMed database {www.ncbi.nlm.nih.gov). 

The C-terminally-truncated human prion protein variant Y145Stop (i.e. HuPrP(23- 
144)) is a valuable model for understanding the fundamental properties of amyloid for- 
mation. To examine the role of AGAAAAGA segment in fibrillization of PrP(23-144), a 
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Figure 1: Identifying the amyloid fibril formation property in PrP(109-136) re- 
gion by the fibril prediction program of l^. 

deletion variant A(l 13-120) PrP(23-144) (in which the palindrome sequence is missing) 
is used Q. The deletion results in an altered amyloid /3-core without affecting amyloido- 
genicity or seeding specificity; this concludes that the core of some amyloids contains "es- 
sential" (nucleation-determining) and "nonessential" regions, with the latter being flexible 
in amino acid sequence space [7 |. 

The amyloid fibrils formed by the polypeptides of PrP(l 13-127), AGAAAAGAVVG- 
GLGG, are taken as the model compound to investigate the biophysical principles gov- 
erning the steric zipper amyloid fibril formation [j4] |5] |8] |9] . The target fibrils adopt the 
structural motif of Class 7 steric zipper, which is formed by stacking of antiparallel /3- 
sheet layers with residue 117 + forming backbone hydrogen bonds to residue 120 — A; 

Computer simulations of amyloid fibril formation by the Syrian hamster prion pro- 
tein (SHaPrP) residues AGAAAAGA, the mouse prion protein (MoPrP) residues VA- 
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Figure 2: The 3D molecular structure of human PrP. (A): the 3D molecular NMR 
structure of human PrP(l 10-136) in dodecylphosphocholine micelles. (B): the 3D 
molecular NMR structure of human PrP( 125-228). 

GAAAAGAV, and their variations GAgG (a longer uninterrupted Alanine stretch flanked 
by Glycine), (AG)4 (a complete disruption of hydrophobic residues), Ag, GAAAGAAA (a 
mimic of A/3(29-36)), Aio, Vio, GAVAAAAVAG (uninterrupted hydrophobic sequence), 
VAVAAAAVAV (less flexible than MoPrP(lll-120)) are studied in [lO]. The first two 
peptides are thought to act as the velcro that holds the parent prion proteins together in 
amyloid structures and can form fibrils themselves (W). 

AGAAAAGA of HuPrP(105-210) is reported being involved into the oligomerization 
and engaged in intra- and/or inter-molecular interactions [ill] . The SHaPrP(109-122) 
peptide near the AGAAAAGA region is observed to form steric zipper fibrils by the data 
of crystal solid-state NMR and molecular dynamics (MD) f51. 

The cellular isoform of the prion protein PrP*^ is located at the cell membrane. Studies 
have shown that exposure of cells to copper (Cu) causes internalisation of PrP*^ in vitro, 
and deletion mutation studies have shown that the palindromic region, amino acids 113- 
120 with the sequence AGAAAAGA is essential for copper-induced internalisation to 
occur |[T2l . Kourie et al (2003) studied the copper modulation of ion channels of PrP(106- 
126) mutant prion peptide fragments and found that the hydrophobic core AGAAAAGA 
is not a Cu^"*" -binding site (but at Met 109 and His 111 ion channels can be formed) [ 13|. 

The AGAAAAGA palindrome in PrP is required not only for the attainment of the 
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PrP^'^ conformation but also to generate a productive PrP^'^-PrP'^ complex that leads to 
the propagation of PrP^'^ llT4l . In contrast to wild type (wt) PrP, PrP lacking the palin- 
drome (PrP A(112-119)) neither converted to PrP^'^ nor generated proteinase K-resistant 
PrP In 1,141. we also know that synthetic peptides corresponding to the so-called "toxic 
peptide" PrP(106-126) segment form fibrils in solution with /3-sheet structure, suggesting 
that PrP(106-126) segment may feature in PrP^'^-PrP*- associations, and larger peptides 
of PrP(107-142) antagonize the in vitro conversion of PrP*^ to the protease-resistant state 
in a cell-free conversion model. 

Zanuy et al (2003) reported that, for AGAAAAGA, the antiparallel strand orientation 
is preferred within the sheets but the parallel orientation is preferred between sheets lITSl 
[161 . AGAAAAGA is one of the most highly amyloidogenic peptides and oligomers of 
AGAAAAGA were found to be stable when the size is 6 to 8 (hexamer to octamer) [il6l . 
Here the AGAAAAGA model of [ 16 1 used for their MD is a homology structure. 

In [III, the following bioinformatics on AGAAAAGA is known. AGAAAAGA dis- 
plays the highest tendency to form amyloid. Peptides containing AGAAAAGA are toxic 
to neurons in culture, whereby the sequence AGAAAAGA was found to be necessary 
but not sufficient for the neurotoxic effect lH]. The synthetic peptides derived from the 
central part of PrP'^(106-141) have an inhibitory effect on this conversion of PrP"= to PK- 
resistant PrP^'^. The presence of residues 119 and 120 (the two last residues within the 
motif AGAAAAGA) seems to be crucial for this inhibitory effect. Mutant PrP molecules 
carrying deletions of amino acids 108-121 or 114-121 are not convertible to PrP^'^; there- 
fore, the central hydrophobic region, spanning all or most of the sequence AGAAAAGA, 
plays an important role in the PrP^'^-PrP*^ conversion process. Wegner et al (2002) as- 
sessed the effect of mutations at and around the AGAAAAGA hydrophobic sequence on 
protease resistance and found that mutations in the central AGAAAAGA hydrophobic 
region lead to immediate alterations in PrP structure and processing 1171 . 

The amyloidogenic and hydrophobic core AGAAAAGA has been implicated in mod- 
ulation of neurotoxicity and the secondary structure of PrP(106-126) ifTSl [T91. which is 
dependent on the formation of aggregated fibril structures regulated by the AGAAAAGA 
core mi. 

Brown (2000) reported that AGAAAAGA blocks the toxicity of PrP(106-126), sug- 
gesting that this sequence is necessary (but insufficient) for the interaction of PrP(106- 
126) with neurons and targeting or use of the AGAAAAGA peptide may represent a 
therapeutic opportunity for controlling prion disease lH] . 

HuPrP(106-126) has been shown to be highly fibrillogenic and toxic to neurons in 
vitro. Jobling et al (1999) found that the AGAAAAGA PrP(l 13-120) hydrophobic core 
sequence is important for PrP(106-126) toxicity probably by influencing its assembly into 
a neurotoxic structure and the hydrophobic sequence may similarly affect aggregation and 
toxicity observed in prion diseases ||T9l . 

Chabry et al (1998) reported that peptides from the central part of the hamster PrP 
106-141 (where residues in the vicinity of positions 106-141 of PrP^'^ and/or PrP*- are crit- 
ically involved in the intermolecular interactions that lead to PrP^'^ formation) could com- 
pletely inhibit the conversion induced by preformed PrP^'^ and the presence of residues 
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119 and 120 from the highly hydrophobic sequence AGAAAAGA (PrP( 113-120)) was 
crucial for an efficient inhibitory effect [20 1. 

Holscher et al (1998) reported that the presence of AGAAAAGA of mouse PrP*^ plays 
an important role in the conversion process of PrP*^ into PrP^'^ and that a deletion mutant 
lacking these codons indeed behaves as a dominant-negative mutant with respect to PrP^'^ 
accumulation ||2l . 

PrP AGAAAAGA is the most highly amyloidogenic peptide and is conserved across 
all species 11211 l22l . Gasset et al (1992) reported there are similarities between the PrP 
sequence AGAAAAGA and that of silkworm fibroin, and the homology between PrP 
sequence AGAAAAGAVVGGLGG and that of spider fibroin ||2T]. Thus, we may say 
that the hydrophobic region of PrP( 109-136) should be a region to form /^-sheets and 
amyloid polymers, instead of a-helices of the Garnier-Robson analysis Il23l . 

2.2 3D Structure of Prion AGAAAAGA Amyloid Fibrils 

Seeing the above survey, we know that prion AGAAAAGA peptide has been reported to 
own an amyloid fibril property (initially described in 1992 by Gasset et al of Prusiner's 
Group). However, there has not been traditional X-ray or NMR experimental structural 
bioinformatics for this octapeptide yet, due to the unstable, noncrystalline and insoluble 
nature of this region, which just falls within the N-terminal unstructured region of prion 
proteins. Studies on atomic-resolution structures of the AGAAAAGA peptide will prove 
useful in future experimental studies on this region, aspects of the structure or the dy- 
namics of this region should play a role in the aggregation process, and knowledge of 
these may be useful for the goals of medicinal chemistry for controlling prion diseases. 
Zhang (201 1) successfully constructed three amyloid fibril models (denoted as Models 
1-3) for the PrP(l 13-120) AGAAAAGA region ||24]. Using all the PDB templates of 
m (with PDB IDs: 20KZ, 20NW, 20LX, 20MQ, 20N9, 20NV, 20NA, lYJO, 20L9, 
20MM, 20NX, 20MP, lYJP), simulation conditions are completely consistent with the 
experimental works of [8] and the simulation methods of package Amber 10 |[25l are: 
the optimization steepest descent (SD) method and conjugate gradient (CG) method to 
relax the Models, then the standard simulated annealing method to make the equilibra- 
tion of Models sufficiently stable and then the SD-CG optimization methods to refine 
the Models. Model 1 and Models 2-3 belong to Class 7 and 1 of [8] respectively, i.e. 
Model 1 is /3-strand antiparallel, face=back, up-up (Figure 3) (where numerical results 
show the agreement with ||4j [l6l), and Models 2-3 are /3-strand parallel, face-to-face, 
up-up (Figures 4-5). In all these models, there is about 5 angstroms between the two 
closest adjacent /3-sheets, maintained by hydrophobic bonds, and about 4.5 angstroms 
between the two closest adjacent /3-strands, which are linked by hydrogen bonds. Illumi- 
nated by PDB templates 3FVA, 3NHC/D, 3NVF/G/H/E, 3MD4/5, 20MP computational 
approaches of global optimization, local search energy minimization (EM), simulated an- 
nealing (SA) and structural bioinformatics etc or introducing novel mathematical formu- 
lations and physical concepts into molecular biology may allow us to obtain a description 
of the protein 3D structure at a submicroscopic level for prion AGAAAAGA amyloid 
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Model 1 (a) Model 1 (b) 



Figure 3: Model 1. (A): blue dashed lines denote the hydrogen bonds between 
the pairs of j3- strands, and the 13 -sheets are maintained by van der Waals contacts 
and hydrophobic packings. (B): the red circle denotes there are very few water 
molecules in this part of the truncated octahedral box of TIP3P waters, the tra- 
jectories of the movement of the amyloid fibril formulate a "cage" made of waters 
and the amyloid fibril flies to the bottom of the "cage", the mouth of the "cage' is 
open - this illuminates to us the amyloid fibril is very hydrophobic. 

fibrils ||2a|26l|27l|28l|29l|30l. 

3 Structural Studies on the PrP(119-131) GAVVG- 
GLGGYMLG 

Some 30 segments from the Alzheimer's amyloid- /3 (A/3) and tau proteins, the PrP prion 
protein, insulin, etc form amyloid-like fibrils, microcrystals that reveal steric zipper struc- 
tures [8|. Harrison et al (2006) reported there are similarities between A/3 and PrP in the 
segment of the three GxxxG repeats (where both A/3 and PrP have the crucial residue 
Methionine located in the middle (GxMxG) of the last repeat) |[3ni32l that controls prion 
formation fS']. Harrison et al (2010) used cell biological approaches of investigating nu- 
merous mutants in this region to reveal the mechanism of prion inhibition, and mutage- 
nesis studies demonstrate that minor alterations to this highly conserved region of PrP*^ 
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Model 2 (a) 



Model 2 (b) 



Figure 4: Model 2. (A): blue dashed lines denote the hydrogen bonds between 
the pairs of j3- strands, and the 13 -sheets are maintained by van der Waals contacts 
and hydrophobic packings. (B): the red half-circles denote there are very few 
water molecules in this parts of the truncated octahedral box of TIP3P waters, 
the trajectories of the movement of the amyloid fibril formulate a "cage" made 
of waters with two open "mouths " circled in red - this illuminates to us that the 
hydrophobic property of amyloid fibrils can form a very interesting movement 
pattern of trajectories. 

drastically affect the ability of cells to uptake and replicate prion infection in both cell and 
animal bioassay 131. This section presents some explanations for the biological experi- 
mental performance of f3 1 from the molecular structural point of view. 

First, we do the alignments of the structured region of mouse, human, dog, rab- 
bit, horse and elk PrP^ (with PDB IDs 1AG2, IQLX, IXYK, 2FJ3, 2KU4, IXYW re- 
spectively) (Figure 6). Figure 6(A) shows the PrP(125-136) 3D structures of human, 
dog, rabbit, horse and elk superposed onto mouse with backbone-atom-RMSD values 
2.579427, 2.228940, 2.745877, 2.532690, 2.877734 angstroms respectively. In Figure 
6(B), the alignment of sequences was generated by the online CLUSTAL 2. 1 program 
at http://www.ehi.ac.uk/Tools/msa/clustalw2/ . where "*" means that the residues or nu- 
cleotides in that column are identical in all sequences in the alignment, ":" means that 
conserved substitutions have been observed, "." means that semi-conserved substitutions 
are observed, the red colored residues are Small (small-i- hydrophobic (incl. aromatic - 
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Model 3 (a) Model 3 (b) 



Figure 5: Model 3. (A): blue dashed lines denote the hydrogen bonds between 
the pairs of j3- strands, and the 13 -sheets are maintained by van der Waals contacts 
and hydrophobic packings. (B): the red half-circles denote there are very few 
water molecules in this parts of the truncated octahedral box of TIP3P waters, 
the trajectories of the movement of the amyloid fibril formulate a "cage" made 
of waters with two open "mouths " circled in red - this illuminates to us that the 
hydrophobic property of amyloid fibrils can form a very interesting movement pat- 
tern of trajectories, and this also illuminates us the property of periodic boundary 
calculations in PMEMD (particle mesh eward molecular dynamics) simulations. 

Y)), the blue colored residues are Acidic, the green colored residues are Hydroxyl + 
sulfhydryl + amine + G, and the gray colored residues are Unusual amino/imino acids 
etc; the residue numbers are from 119 to 231 (for rabbit the numbers are 1 less than 
others); the last column of numbers denotes the number of residues accounted from 1; 
the PrP 125-136 LGGYMLGSAMSR, /3-strand 1 (/31), a-hehx 1 (al), /3-strand 2 (/32), 
a-helix 2 {al) and a-helix 3 (a3) of a PrP structure (Figure 2(B)) were underline de- 
noted. In Figure 6, we can see the PrP( 125-136) LGGYMLGSAMSR is highly conserved 
among all species, strongly suggesting it has functional and evolutionary significance. To 
understand the functions of proteins at a molecular level, in protein structures, the non- 
covalent interactions such as hydrogen bonding (HB), ionic interactions (SB), van der 
Waals forces (vdW), and hydrophobic packing (HP) are driving the proteins to be able to 
perform their biological functions {http://en.wikipedia.org/wiki/Protein_stnicture). Thus, 
in the below, we will investigate the HBs, SBs, vdWs, HPs in the structures and their 
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CLUSTAL 2.1 multiple sequence alignment 



ARG136 




PrP( 125-136) 



(a) 



ZFJ3_rabbit.pd]o 

IXYWelk.pdb 

IXYKdog.pdb 

2KU4_horEe.pdb 

lAG2_niouse.pdb 

IQLX humsn.pdb 



2FJ3_rabbit.pdb 
lXYW_ellc.pdlj 
lXYK_dQg.pdb 
2KU4_horEe.pdb 
lAG2_mQuse .pdb 
lQLX_human . pdb 



(b) 



l.GGTMI.GSAWSR PLIHFGNDYEDRYYREWMYRYPNQVYYRPVDaYSMQNSFVHD 54 

PNQVYYRPVDQYNHQNTFVHD 58 
PDQVYYRPVDQYSNQWHFVRD 58 
PHQVYYRPVSEY3NQKNFVHD 60 
■PNQVYYRPVDQYSMQHNFVHD 55 
PHQVYYRPMDEY3HQNNFVHD 54 



--WGG LeGYMLGSAMSR PLIlIFGMDYEDRYYREHMYRYP 
■ -WGG LGGYMI.G5AM5R PLIHFGMDYEDRYYREHMYRYP 
GSWGG LEGYMI.G3AMS RPLIEFGMDYEDRYYREHMYRYP 
GLGGYMLGSAMSRPMIHFGMDWEDRYYREHMYRYP 



LSGYMLGSaMSR PIIHFGSDYEDRYYREHMHRY 

9 01234567 8 9012 34 5678 90123 4567 89 0123 4 567 
119 130 140 150 

pi al 



901234567 8 9012 345678 
160 170 
fi2 



CVNITVKQHTVTTTTKGENFTETDIKIMERWEQMCITQYQQESQAAYQRA-- 105 
CVNITVKQHTVTTTTKGENFTETDIKMMERWEQMCITQYQRESEAYYQRGA3 111 
CVNITVKQHTVTTTTKGENFTETDMKIMERWEQMCVTQYQKE3EAYYQRGA3 111 
CVNITVKQHTVTTTTKGENFTETDVKIMERVVEQMCITQTQKEYEAFQQRGA3 113 
CVNITIKQHTVTTTTKGEMFTETDVKMMERVVEQMCVTQYQKESQAYY 103 
CYNITIKQHTVTTTTKGENFTETDVKMMERWEQMCITQYERESQAYYQR 104 



(-1) 



************* 



901234567890 123456789 0123456789012345678901234567 8 901 
180 190 200 210 220 231 



Figure 6: The alignments of the structured region of mouse, human, dog, rab- 
bit, horse and elk proteins with PDB IDs 1AG2, IQLX, IXYK, 2FJ3, 2KU4, 
IXYW respectively. (A) shows the 3D structural alignments of PrP^ (125-1 36) 
(whose residues were underlined in (B)). (B) shows the ID sequence alignments 
of PrpC( 119-231). 

structural dynamics of LGGYMLGS AMSR of human and rabbit, GLGGYMLGSAMSR 
of mouse, VVGGLGGYMLGSAMSR of elk and dog, GS VVGGLGGYMLGS AMSR of 
horse (where S120 is special for horse, instead of A120) and structural connections with 
other residues/loops/sheets/helices in the C-terminal. Mutations will destroy these non- 
covalent interactions that well maintain the structure so the function of the prion protein. 
This will give clear explanations to the mutants in the Gly cine-rich region of E l33l . 
which affect the uptake of prion infectivity very much. 

In 2010, horses were reported to be resistant to prion diseases fi^. First, we analyze 
the role of GSVVGGLGGYMLGSAMSR(PrP(l 19-136)) in horse PrP'^( 119-231) (PDB 
entry 2KU4). In our study |[35l . the MD simulation conditions are at 350 K in explicit sol- 
vent under neutral and low pH environments, heatings are using the Langevin thermostat 
algorithm in constant NVT ensembles and the equilibrations and productions are using 
Langevin thermostat algorithm in constant NPT ensembles. Seeing the Tables I-III of the 
Supplementary Material of El, we know the following HBs of GSVVGGLGGYML- 
GSAMSR(PrP( 119-136)): 

ARG136-TYRI57 (linking /31-to-al loop with al-to-/32 loop, 68.27%, 71.04%), 

MET134-ASN159 (linking /31-to-al loop with al-to-/32 loop, 29.83%, 21.47%), 

GLY131-GLN160 (linking /31 with /32, 26.67%, 37.03%), 

SER132-GLN217 (linking /31 with a3, 29.70%, 12.90%), 

ARGI36-PR0158 (linking /31-to-al loop with al-to-/32 loop, 9.53%, 6.01%), 
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ARG136-TYR157 (linking /31-to-al loop with a\-to-/32 loop, 68.27%, 71.04% ), 
GLY126-ARG164 (linking bend before /31 with /32-to-a2 loop, 33.31%^), 
TYR128-ASP178 (hnking coil before /31 with a2, 11.81%^), 
SER120-LEU125 (in the peptide, 9.72%^), 
LEU125-10TYE128 (in the peptide, 5.58%2), 
GLY119-4VAL122 (in the peptide, 6.51%2), 

GLY127-ARG164 (link bend before /3\ with ^2-to-a2 loop, 5.87%^), 

where the first percentage is for seedl (%^) and the second percentage is for seed2 (%^) 
and the two seeds mean two different initial velocities of 30 ns MD, and the following 
HPs of GSVVGGLGGYMLGSAMSR (PrP(l 19-136)): 

In PrP(l 19-136) of 2LBG.pdb, between the two adjacent residues there are always 
occupied by HPs with rate of 100%, except for between GLY123 and GLY124, and 
between GLY126 and GLY127, 

Among the residues in PrP( 119-136), there are HPs in 

TYR128-LEU130 (where LEU130 is a residue in /31), 

VAL121-GLY119, GLY123, TYR128, MET129, LEU130 (where MET129 
and LEU130arein /31), 

VAL122-SER120, GLY124, LEU125, TYR128, 
GLY123-SER120, LEU125, 
TYR128-GLY126, 

MET129-GLY131 (where both MET129 and GLY131 are in /31), 
MET134-SER132, ARG136 (where SER132 is in /31), 

The HPs between a residue in PrP(l 19-136) and a residue in PrP(137-231): 

SER135-PR0137 (where SERBS, PRO 137 is in /31-to-al loop), 

ARG136-MET154, TYR157 (where ARG136 is in /31-to-al loop, MET154 
and TYR157 are in al-to-/32 loop), 

ASN159-SER135, ALA133, ARG136, MET134 (where ASN159 is in al- 
to-/32 loop), 

TYR162-TYR128, MET129, LEU130 (where MET129, LEU130 are in ^1, 
and TYR162 is in /32, there is no HB between TYR162 and LEU130), 

GLN217-GLY131 (where is in /31, GLN217 is in a3). 
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ARG164-GLY127, MET129, TYR128 (where MET129 is in (31, ARG164 
is in /32-to-a2 loop), 

TYR163-TYR128, LEU130, MET129 (where MET129, LEU130 are in /31, 
TYR163 is in /32), 

GLN162-ALA133, LEU130, MET134, GLY131, SER132 (where ALA133, 
LEU130, GLY131, SER132 are in /31, GLN162 is in /32), 

VAL161-GLY131, LEU130 (where GLY131, LEU130 are in VAL161 is 

in /32), 

MET213-MET134 (MET213 is in a3, MET134 is in the /31-to-al loop), 

with > 50% occupancy rate over the long MD ti^ajectory of 30 ns. The mutations listed in 
|[3l will break these HBs and HPs to lose the inhibition to prion diseases of horses. 

Rabbits are also resistant to infection from prion diseases from some species ||36] l37l 
[381 [391 SOI and the outbreak of "mad rabbit disease" is unlikely EH Sll- Nisbet et al 
(2010) reported there is 87% sequence homology between mouse PrP and rabbit PrP, ap- 
proximately 9 of the 33 (i.e. 1/3) of the difference in the region DGRRSSSTV of mouse 
and QRAAGVL of rabbit, and residues surrounding the glycosylphosphatidylinositol an- 
chor attachment site of PrP modulate prion infection ||431 . Under the MD simulation 
conditions of 450 K in explicit solvent, neutral and low pH environments, with heatings 
using the Langevin thermostat algorithm in constant NVT ensembles and the equilibra- 
tions and productions using Langevin thermostat algorithm in constant NPT ensembles, 
we have found rabbit PrP has the C-terminal residue R227 forming a HB/SB network 
with inner residues and the beginning residues of rabbit homology strcture PrP( 120-229) 
(6EPA.pdb) and NMR structure PrP(124-228) (2FJ3.pdb) but mouse PrP has no such an 
Arginine residue at the end of C-terminal owning this property Il44ll45]| . 

Dogs were reported in 2008 to be resistant to infection from prion diseases from other 
species f46l. In our recent work [47 1, MD studies were done under 450 K in explicit sol- 
vent, neutral and low pH environments, with heatings using the Langevin thermostat algo- 
rithm in constant NVT ensembles and the equilibrations and productions using Langevin 
thermostat algorithm in constant NPT ensembles. In the hydrophobic region 121-136 of 
dog PrP, MD studies find that there are strong HBs S132-Y163-Y128 linking the two 
antiparallel /3-strands, and strong HPs M134-A133, M129-L130, V121-V122 residing 
in the core of the hydrophobic region 121-136 PTI . We should note that residue 129 for 
dog PrP is LI29 (http ://www. uniprot.org/unrprot/Q4650 1 ) but for many others is M129. 

For the hydrophobic region PrP( 109-136), we should also notice the following points: 

(1) the segment GYMLGS or GYVLGS of HuPrP(127-132) can form Class 8 an- 
tiparallel amyloid fibrils (3NHC.pdb and 3NHD.pdb) |[8ll48ll. 

(2) the mutants A117V and M129V cause GSS EH [50l [SB and CJD |[52l l53l 

respectively. 
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(3) the /3-sheet core of PrP^'^ consists of three layers of /3-strands El(116-119), 
E2(129-132) and E3(160-164) ||54l, where El and E2 are in the hydrophobic re- 
gion (109-136), 

(4) HI 1 1 is a residue of copper binding sites Il55ll56]l57ll58) . 

(5) Y128 is in the center of HB-and-SB -network of HB(Y128-D178), SB(D178- 
R164), HB(Y128-R164), HB(Y128-H177), HB(H177-N154) f59ll60ll6ni52ll. 

(6) A 133 and SI 32 have HBs with R220 and a water binding site with G131 162 
l64l . in PrP(l 13-132) the hydrophobic cluster with vdWs rendering of atoms in 
residues 113-127 interacts with the first /3-strand (see Figure 2(B) of [^5l), M129 
makes interactions with the side chain of V122 and pulls the N-terminus into the 
/3-sheet [611 and M129 is very close to Y163 ED, 

(7) conservation of the Gly-rich region PrP(l 19-131) is required for uptake of prion 
infectivity |l3]|66j|67l is shown by the physical or chemical properties of numerous 
mutations in the PrP(l 19-131) Gly-rich region ||3], 

(8) one 0-linked sugar at Serl35 can affect the coil-to- /3 structural transition of the 
prion peptide, but at Serl32 the effect is opposite ll68l . etc., and 

(9) the NMR structure of HuPrP( 1 10-136) in dodecylphosphocholine micelles was 
known (2LBG.pdb) ||69ll and we do MD simulations for it as follows. 

We used the ff03 force field of the AMBER 11 package in a neutral pH en- 
vironment. The systems were surrounded with a 12 angstroms' layer of TIP3PB0X 
water molecules and with 2 CI- ions added using the XLEaP module of AMBER 11. 
The templates used are the 2LBG.pdb from the Protein Data Bank and its 17 mutants 
at G114V, A117V, G119A, G119L, G119P, A120P, G123A, G123P, G124A, L125A, 
G126A, G127A, G127L, M129V, G131A, G131L, G131P, which are got by the mu- 
tate module of the free package Swiss-Pdb Viewer Version 4.1.0 (http://spdbv.vital-it.ch7]l. 
These 18 models were firstly optimized by SD method and then CG method for 2 stages. 
Minimization Stage 1 is holding the solute fixed with a force constant of 500.0 kcal mol"' 
angstrom"^ for 500 steps of SD minimization followed by 500 steps of CG minimization. 
Minimization Stage 2 is minimizing the entire system for 1000 steps of SD minimization 
followed by 1500 steps of CG minimization. The minimized models were checked by 
Swiss-Pdb Viewer and found there is not any bad contact that makes amino acids clashed. 
Then, the solvated proteins were quickly heated from K to 310 K linearly for 300 ps 
and then systems were kept at 310 K for 700 ps. The systems were in constant NVT 
ensembles using Langevin thermostat algorithm with weak restraints (a force constant of 
10.0 kcal mol"' angstrom'^) on the solvated proteins. The SHAKE and PMEMD algo- 
rithms with nonbonded cutoffs of 12 angstroms were used. Next, the systems were done 
MD simulations in constant NPT ensembles (with constant pressure 1 atm and constant 
temperature 310 K) under the Langevin thermostat for 4 ns and the PRESS, VOLUME 
(DENSITY) and RMSDs were sufficiently stable for each of the 18 models (Figure 7). A 
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step size of 2 fs was used for all the MD simulations, the structures were saved to file every 
1000 steps and the Metropolis criterion was used. These MD simulation conditions are 
completely consistent with the experimental work of NMR structure of HuPrP(l 10-136) 
(2LBG.pdb). 




Figure 7: Variations of the potential energy (EPTOT), kinetic energy (EKTOT), 
total energy (ETOT), temperature, pressure, volume, density, and RMSDs during 
the 5 ns ofMD simulations for all the models. In the RMSD graph, different colors 
stand for HuPrP( 110-136) and its 1 7 mutants. 

We picked out the snapshots at 0, 1, 2, 3, 4 and 5 ns for each model (Figure 8) and 
found the a-helix structure of each model has been unfolding and longer MD simulations 
might make the a-helix structure unfolded completely (Figures 8~9). This shows that 
the structure in PrP(110-136) region is structurally unstable and might be critical to the 
conversion from the predominantly a-helical PrP*-^ into the rich /3-sheet PrP^'^. 

The preliminary findings to explain the above performance during the 5 ns of MD 
simulations may be described as follows. At the beginning of the simulations, the 3D 
structures of all the 13 models are a-helices (Figures 8~9). with rich HBs. In an ideal a- 
helix, there are 3.6 residues per complete rotating so a rotation of 100 degrees per residue. 
Thus, in all the 13 models, there are about 7 turns (Figure 8) in each a-helix. A a-helix is 
maintained by HBs. Clearly, the disappearance of HBs is a reason for the unfolding of all 
the 13 models. As the progress of MD simulations, most of them have disappeared, except 
for the three main HBs GLY/ALA/LEU/PR0131-SER135-ARG136, TYR128-SER132, 
which have different occupied rate for each model (Table 1) and are just in the N-terminal 
of the structural region of HuPrP(125-228) (IQLX.pdb). This might show that PrP(110- 
124) region of PrP(l 10-136) (2LBG.pdb) is very unstable. We also found that there is 
a SB between HISlll-LYSllO with the occupied rate 100% for the 13 models. Spe- 
cially for the mutants G127L, M129V, G131A, and G131L, there is another SB between 
HIS111-136ARG linking the head and tail of HuPrP(l 10-136). Seeing the snapshots 
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Figure 8: The respective snapshots at ns, 1 ns, 2 ns, 3 ns, 4 ns and 5 ns of the 
MD simulations for the 18 models. The dashed lines denote hydrogen bonds. 

of 3, 4 and 5 ns of mutant G127L in Figure 8, 4 and 5 ns of mutant M129V in Figure 
8, 3, 4 and 5 ns of mutant G131A in Figure 8, and 3, 4 and 5 ns of mutant G131L in 
Figure 8, we may know the SB HIS111-136ARG makes these snapshots looking like a 
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Figure 9: Variations of the secondary structures of the 18 models. H is the a-helix, 
B is the residue in isolated (i-bridge, E is the extended strand, G is the 3-helix or 
3/10 helix, I is the it -helix, T is the hydrogen bonded turn and S is the bend as 
described in the DSSP program IITTJ . 

"hairpin". HISlll is a very important residue in PrP(l 10-136). Along with the unfold- 
ing of a-helical structure of all these 13 models, we found many HPs disappeared except 
for some fundamental ones MET134-ALA133, LEU 1 30-1 29MET, VAL122-VAL121- 
ALA120, ALA118-ALA117-ALA116-ALA115, and ALA113-MET112 with the oc- 
cupied rate 100%, where ALA118-ALA117-ALA116-ALA115 are in the core of the 
palindrome AGAAAAGA and this might imply to us the hydrophobic core is very hard 
to break and this palindrome really has enormous potential to be amyloid fibrils. 

4 Concluding Remarks on PrP(109-136) 

To really reveal the secrets of prion diseases is very hard. For us it is a long shot but 
certainly worth pursuing. It was reported that the hydrophobic region PrP( 109-1 36) con- 
trols the formation into diseased prions: the AGAAAAGA palindrome and Glycine-xxx- 
Glycine repeats (both being the inhibitor of prion diseases) are just in this region. This 
paper gives some investigations and explanations on the PrP( 109-136) region in view of 
its 3D structures and molecular dynamics studies. The structural bioinformatics presented 
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Table 1: Occupied rate (%) of the hydrogen bonds for each model 



HuPrPtl 10-136) 


G114V 


A117V 


G119A 


131@0-135@OG.HG (25.48) 
125@0-135@OG.HG (17.68) 

1 17(3)0—1 "^fifffiJ^JT-TI HHI 1 (1 '^')\ 

116@0-136eNHl.HH12 (5.76) 


131@0-135®OG.HG (26.28) 
135@0-136@NH1JIH11 (21.72) 

132@0-135@OG.HG(7.72) 


131@0-135®OG.HG (40.52) 
128@0-132@OGJIG (8.88) 

132®0-135@OG.HG (7.36) 


131@0-135®OG.HG (49.24) 
135@0-136®NH1JIH11 (14.6) 


G119L 


G119P 


A120P 


G123A 


135@0-136@NH1JIH11 (23.96) 
131@0-135@OGJIG (22.8) 
132@0-135@OG.HG (15.08) 


131@0-135@OG.HG (31.24) 
135@0-136@NH1.HH11 (24.4) 
128@0-132@OGJIG (10.08) 
130@O-135@OG.HG (7.2) 


135@0-136@OG.HG (28.08) 
132@0-135®OG.HG(8.6) 
135@0-136@NH1JIH11 (7.72) 
112®0-111®ND1.HD1 (5.6) 
128®0-132@OGJIG (5.04) 


131@0-135@OG.HG (28.08) 
135®0-136®NH1.HH11 (17.76) 
128®0-132®OGJIG (6.16) 


G123P 


G124A 


L125A 


G126A 


131@0-135@OG.HG (25.52) 
13S@0-136®NH1.HH11 (21.92) 


13ie>0-13S@OG.HG (41.4) 
13S@'0-136<aiNHl.HHll (21.4) 

128®0-132@OGJIG (5.56) 
134@0-136@NH1.HH11 (5.44) 


131@'0-135@OG.HG (36.24) 
135@0-136eNHl.HHll (17.92) 
132®0-135®OG.HG (8.48) 
128®0-132@OG.HG (5.6) 


13ieO-135eOG.HG (46.36) 
135eO-136eNHl.HHll (14.12) 
128eO-132@OG.HG (7.04) 
133eO-136eNE.HE (5.6) 
132<a>0-135(a>OG.HG (5.08) 


G127A 


G127L 


M129V 


G131A 


131@0-135@OG.HG (26.2) 
135@0-136®NH1JIH11 (14.88) 
112@0-111@ND1.HD1 (13.08) 


13S@0-136@NH1.HH11 (23.6) 
131@0-135®OG.HG (21.2) 
133®0-135®OG.HG (8.04) 


131@0-135®OG.HG (34.4) 
135@0-136®NH1JIH11 (13.2) 
128@0-132®OG.HG (8.92) 


13ieO-13SeOG.HG (47.76) 
128@0-132eOG.HG (17.08) 

135@0-136®NH1.HH11 (13.48) 
132@0-136eNHl.HHll (7.72) 

134®0-111HIS®ND1.HD1 (6.2) 




G131L 


G131P 






131@0-135®OG.HG (70.4) 
135@0-136@NH1JIH11 (18.4) 
128@0-132®OG.HG (7.88) 


131@0-135@OGflG (32.48) 
135@0-13(>@NH1.HH11 (16.72) 
128@0-132®OG.HG (15.2) 





in tills paper can be acted as a reference in 3D images for laboratory experimental works. 
This presents some clue or hints for the author to study prion proteins and prions. 
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